Island nucleation in the presence of step edge barriers: 
Theory and applications 
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O , We develop a theory of nucleation on top of two-dimensional islands bordered by steps with 

• an additional energy barrier AEs for descending atoms. The theory is based on the concept of 

the residence time of an adatom on the island, and yields an expression for the nucleation rate 
which becomes exact in the limit of strong step edge barriers. This expression differs qualitatively 
■ and quantitatively from that obtained using the conventional rate equation approach to nucleation 

[J. Tersoff et ai, Phys. Rev. Lett. 72, 266 (1994)]. We argue that rate equation theory fails 
because nucleation is dominated by the rare instances when two atoms are present on the island 
' simultaneously. The theory is applied to two distinct problems: The onset of second layer nucleation 

in submonolayer growth, and the distribution of the sizes of top terraces of multilayer mounds under 
conditions of strong step edge barriers. Application to homoepitaxial growth on Pt(lll) yields the 
estimate AEs > 0.33 eV for the additional energy barrier at CO-decorated steps. 
PACS: 68.55.-a, 81.15.-z, 68.10.Jy, 05.40.-a 



I. INTRODUCTION 



On many metal surfaces an adatom descending across a step edge, encounters an additional energy barrier AEs 
compared to the diffusion barrier Ed on an atomically flat terraced. This step edge barrier controls, the rate of 
interlayer mass transport and therefore has a decisive influence on the morphology of multilawij-SlmsBcl. 

A number of papers have addressed the question of how to determine AEs experimentallyullil. Most approaches 
start from the observation that the step edge barrier increases the residence time of adatoms on top of an island, 
and therefore promotes the nucleation of the next layer. Tersoff, Denier van derXJon and TrompB (TDT) presented 
a quantitative analysis of this effect, which is based on the conventional estimate!! 2 ] 
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y—{ • for the nucleation rate ui in terms of the in-layer hopping rate v, the adatom density n and the size i* of the largest 
0^ | unstable cluster. Equation (Q) arises from rate equations for spatially averaged island and adatom densities, and its 
0^ i applicability to the confined geometry on top of an island is not obvious. Other treatments invoke the concept of a 
critical adatom density for nucleation, which is defined as the density at which nucleation would take place on the 
unbounded terraceo Ll. Ag ain, given the rather different conditions on top of an island, this approach seems hard to 



justifyd (see Section III C ) 



In this paper we present a detailed microsopic analysis of nucleation in the presence of strong step edge barriers, 
which takes into account the large fluctuations of the adatom population on top of the island. Our expression for 



the nucleation rate, to be derived in Section II A, differs qualitatively and quantitatively from that obtained by 
TDT. The rate equation treatment significantly overestimates the nucleation probability, and therefore the analysis 
of experimental data based on the TDT theory generally gives values for AEs which are too small. To illustrate our 



point, we apply our approach to two different experimental situations. In Section III we determine the critical radius 



for second layer nucleationjn submonolayer deposition, and use it to reanalyse the experiments of Bromann et alia for 



x. 

Ag(lll), and of Kalff et alJ13 for CO-contaminated Pt(lll). In Section IV -Hie. develop a simple analytic theory for the 
size of the top terrace of pyramidal multilayer mounds ("wedding cakes^otj), and apply it to the case of Pt(lll)Ej. 
For CO-decorated steps on Pt(lll) the two approaches yield mutually consistent estimates for AEs- Some of the 



results presented in Section III have been independently obtained by Rottler and MaassliZl, who also provide numerical 



confirmation by Monte Carlo simulations for i* = 1. 



II. NUCLEATION RATE ON TOP OF AN ISLAND 
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A. Irreversible aggregation with strong barriers 



We consider a compact two-dimensional island of area A and perimeter L. Both A and L are dimensionless quantities 
measured in terms of the number of lattice sites on the island (A) and the number of edge sites (L), respectively. For 
large L they are related through 

A = aL 2 (2) 

where the coefficient a depends on the island shape; for example, a = 1/16 for square islands, a — 1/12 for regular 
hexagons on a triangular lattice and a = l/Air for a circular island in a continuum approximation. Atoms are deposited 
at rate F, diffuse on the island (and on the underlying terrace) at rate v = vq exp[— En/hsT], and descend from edge 
sites at rate v' = v'^ cxp[— Es/ksT]. The additional energy barrier at the-sipa edge is then AE$ = E$ — Ed. Since 
the actual energy landscape on real islands may well be more complicated! 1 3cH, v 1 should be regarded as an effective 
interlayer transport rate. r— , 
The problem involves three time scales: The traversal timeE! 

Ttr « A/v (3) 

required for an atom to visit all sites of the island, the mean time interval 

At = {FA)' 1 (4) 

between subsequent deposition events on the island, and the mean residence time t that a single adatom spends on 
the island. The latter is quite generally given in terms of the mean adatom density n and the flux F as 

r = n/F. (5) 

To arrive at this relation, consider the deposition of noninteracting atoms onto the island during some long time 
interval T. A total of T/ At atoms arrive on the island, each of which spends an average time r. The total time spent 
by all the atoms is Tr/At, and hence the mean number of atoms simultaneously present on the island is r/ At = nA, 
from which (||) follows. 

Here-we focus on the case of strong step edge barriers (case 2 of TDT), in the sense that the Ehrlich-Schwoebel 
lengtbB 

^ ES = {y/y' - 1) w v/v> (6) 

is much larger than the island size. In this limit the adatom density on the island becomes spatially uniform!, and its 
value is determined simply by the balance between the number of atoms FA deposited on the island, and the number 
of atoms Lnv' leaving it per unit timeoa. Thus h = FA/Lv' = aFL/v' and 

T = — • 00 

v' 

It follows that t/txt ~ ^es/L, and therefore the strong barrier condition is equivalent to r/r tr ^> 1. 

In this section we further assume that the critical island size takes its minimal value i* = 1, i.e. dimers are stable 
and immobile. Then, as soon as two adatoms are present on the island, nucleation occurs after a time of the order 
of rtr t; on the scale of the residence time nucleation is an instantaneous event. We conclude, therefore, that 
the nucleation probability p nuc per deposited adatom is equal to the probability P2 that two adatoms are present 
simultaneously on the island. To compute P2, suppose the first adatom has arrived on the island at time t = 0, and 
denote by t\ and ti the departure time of the first adatom and the arrival time of the second one. Since deposition is 
a Poisson process, ti is an exponential random variable with mean At. In the strong barrier limit also the distribution 
of residence times t\ decays exponentially^ for times large compared to r tr . Thus we have 

Pnuc = P2 = Prob[h >t 2 ] = — dtt e-^' T / dt 2 e- t2 ' At = — — . (8) 

T~At J J T + At 

If At <C r we find p 2 ~ 1, i.e. the first adatom deposited on the island already forms a nucleus. In this limit the 
nucleation rate becomes independent of v' , which means that the step edge barrier i s, at l east as far as second layer 



nucleation is concerned, infinite. We will return to this case at the end of Section III A. In the following we will 
mostly be interested in the opposite limit At ^> t, where p 2 ~ t/ At <§; 1. A general expression for the probability 
Pn+i to find n + 1 adatoms simultaneously on the island is derived in Appendix [A|. 
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The expression for p nuc can also be found by using an argument first developed in one dimension by Elkinani 
and VillainEfl (see Ref. 25 for the two-dimensional case). According to this approach, the nucleation probability per 
deposited atom is simply the number of distinct sites visited by each adatom (equal to A) times the probability that 
a site is occupied (given by the adatom density n). Together with (JsJ) this implies p nuc = AFt = r/At as above. 

To obtain the nucleation rate to, defined as the number of nucleation events taking place on the island per unit 
time, we multiply the nucleation probability p nuc by the number of atoms arriving on the island per unit time, and 
obtain, in the relevant case At t, 



, t a 3 F 2 L 5 

U = At = ~v' ' < 9 > 



This is to be compared to the expression 



o?F 2 L^v . . 

" = ( 10 ) 

derived from the approach of TDT (7 denotes a capture number of order unity). Equation (|l0| ) exceeds (|^) by a 
factor (v/v')^ 1 = £es/L ^> 1. To explain this discrepancy we note that the mean number of adatoms on the island 
is much less than unity in the regime of interest. Most of the time the island is empty, and sometimes it is occupied 
by a single adatom. For example, consider a Pt island on Pt(lll) without a second laye r nuc leus, but at a size, at 



which about half of the island population has already formed such a nucleus (see Section [II E ). Even for this island 
the probability to find an atom on top is only about 10 -2 , as is the probability for an atom to become deposited on 
an island already populated by an atom. For this reason we refer to our approach as the lonely adatom model (LAM) . 
The strong fluctuations in the occupancy of the island imply that the replacement of the actual adatom density by 
its time average n, which is inherent in ([!]), cannot be justified, and must be replaced by the more detailed statistical 
analysis provided above. 

B. Intermediate and weak barriers 

To some extent the considerations of the preceding section can be generalized to the case when the strong barrier 
condition £es 3> L is no longer satisfied. The mean residence time is still given by (||), and the mean adatom density 
n can be obtained by solving the stationary diffusion equation with the appropriate boundary conditionsEoE£l. For 
the case of a circular island of radius r — L/2tt this yields 

'-5(^7)- <"> 

However now the probability P2 provides only an upper bound on the nucleation probability p nuc , since one of the 
two adatoms may escape from the island before the two meet. Also in the evaluation of P2 in eq.(||) the probability 
distribution P lcs {ti) of residence times has to be used, which is a complicated non-exponential functionO. The 
calculation simplifies by noting that for weak barriers surely r <C At; then the exponential distribution of interarrival 
times can be expanded and one obtains 



f°° 1 f°° T 

P2= J dh P re8 (ti)(l ~ e- l ^ At ) ~^ t J dh hP Ies (h) = 

We conclude that an upper bound on the nucleation rate for a circular island is given by the expression 

■K 2 F 2 r 4 / o vr\ 



(12) 



2v V v / 

If the barriers are weak, r 3> ^eSj we find cj> = ir 2 F 2 r 6 /2u, which agrees (up to a numerical prefactor) with Eq. (4a) of 
TDT. This implies that the rate equation ansatz (^ yields the correct result only when the barriers can be neglected. 

C. Reversible aggregation 

We return to the case of strong step edge barriers and analyze reversible aggregation with stable trimers and 
unstable dimers (i* = 2). The dimer dissociation time Tdi S introduces a further time scale into the problem. Suppose 
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that two adatoms are present on the island; this is true with probability p2 = r/At. Let us assume for simplicity 
that one atom is fixed and the other diffusing. If the two adatoms stay on the island a time t, this means that the 
diffusing atom has failed to descend a number of times of order j = i/(r tr + Td; s ), which is true with a probability 
of the order of order exp(— jr tr /r). Therefore the effective residence time r' of the adatom pair can be estimated as 
t' ~ T(r tr + Tdi S )/7t r . In the limit r tr 3> Tdi s , t' ~ r, and in the limit T tr <C Tdi S , t' ~ T(Tdi S /T tr ). The lifetime of the 
pair is increased compared to the single adatom lifetime if r tr <C Tdi S . 

Consider first the case Tdis Ttr, when the total residence time of the two adatoms is still of the order of t, and a 
dimer is present a fraction q = Tdi s /T tr of that time. The probability for a third atom to deposit while the two atoms 
are on the island is r/At as before. The third atom traverses the island of the order of m — t/t^ times, and each 
time the probability that it encounters a dimer is q. We therefore have to distinguish the cases mq <C 1 and mq 3> 1. 
In the first case the probability that a stable trimer forms is mq ~ TTd; s /T t 2 r , and altogether the nucleation probability 
per atom is of the order of 

/ T \ 2 TT dis F 2 L 3 lS 2 T d 



Pnuc ~{~Ai) ~ 2 ^73 for r di S < T tr/ T (regime I). (14) 

On the other hand, if mq ^> 1 the third atom is certain to encounter a dimer, and 

/ t \ 2 F 2 L 6 

Pnuc ~ (^J for t 2 /t < T dis < Ttr (regime II). (15) 

In fact a somewhat more precise statement can be made. In regime II the nucleation probability is equal to the 
probability p 3 of finding three adatoms simultaneously on the island. The calculation in Appendix |A| shows that, for 
t < At, p 3 = (1/2)(t/A£) 2 , and therefore the prefactor in @ is a 4 /2. 

Next consider the case Td; s 3> T tr , when the effective residence time of the first two adatoms is t' ~ T(Tdi S /T tr ) 3> t. 
If a third atom is deposited onto the island it will find the dimer with a probability close to unity. A third atom will 
deposit during time t' with probability t' /At if t' < At and with probability one if t' > At. This implies two further 
scaling regimes: 

Pnuc ~ ~rr ■ ~t~7 ~ F L ,n Td ' 8 for T tr < T dis < At ■ (t^/t) (regime III). (16) 
At At i/ z 



t FL 3 

Pnuc ~ ^ ~ -j-r- for Tdis > At • (T tI /T) (regime IV). (17) 

In the last regime the dimer lifetime is so long that the situation effectively reduces to the case i* = 1. 

The smallest value that Td; s can take is the inverse in-layer hopping frequency \jv. Then Tdi s /Tt r ~ L~ 2 <C 1 
always, but depending on the strength of the step edge barrier both regimes I and II can be realized, corresponding 
to L 3 ^> £es (regime I) and L 3 <C ^es (regime II), respectively. Setting Tdi S = \ jv in regime I (Eq. fll4j)) yields 



F 2 L 3 v 

Pnuc ~ tt— (regime I, T d i s = l/v) (18) 

h ' ' 

which agrees with the expression derived by TDT for i* — 2. This indicates, as suggested by Rottler and Maass0, 
that the rate equation approach describes a situation in which the probability for an encounter between the atoms 
forming the stable nucleus (the quantity mq in our derivation) is small compared to unity. 

The complexity of the case i* = 2 illustrates the difficulty of extending the present analysis to larger values of i*. 
The generalization seems straightforward only in regime II, where the barriers are sufficiently strong to guarantee the 
formation of a stable nucleus once i* + 1 atoms are on the island simultaneously (but still so weak that t <C At). 
Then p nuc = Pi*+i, and using Eq. ( |A2| ) of the Appendix we obtain 

1 / t y* 1 fa 2 L 3 F\ 1 ' 
Kmc = — { [-^) = I — — J (regime II,generaH*). (19) 



III. SECOND LAYER NUCLEATION 
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A. The critical island size 



In this section we consider a population of submonolayer islands and ask for the fraction / of islands on which a 
second layer has nucleated. Once the nucleation rate ui{L) per island is known as a function of island size L, this 
fraction is given byo 



/ = l-exp[- j dtLo{L(i)% 



(20) 



where the integration extends over the deposition time, and the functional dependence L(t) depends on the growth 
situation. For example, if first layer islands of density N nucleate at time t = , and no further nucleation occurs at 
later times, then the island size increases with coverage 9 = Ft as L — ^/9/aN. Assuming a size dependence of the 
nucleation rate of the general form 



Lo(L) = FflL k , 
where Q is a constant independent of L, one finds 

/ = 1 - cxp[-(L/L c ) fc + 2 ] 

where the critical island size L c is given by 

_ (k + 2 V /(fc+2) 
Lc ~\2aNn) 

In particular, using our expression (m for i* — 1 we obtain 



7 v 1 \ 
2a 4 ' FN J 



1/7 



(21) 
(22) 

(23) 
(24) 



This agrees with the scaling law derived by Rottler a nd M aassEZI and specifies the prefactor. 

The different regimes for i* = 2 discussed in Section II C can be treated in the same way; the results are summarized 
in Appendix ||. Inspection shows that our regime II corresponds to regime II of Rottler and MaassO, while our regime 
I corresponds to their regime III (they assume that the dimers are maximally unstable in the sense that r^is = ^/v). 

A critical island size for second layer nucleation can also be defined when the barriers are so strong (or the flux so 
large) that the condition r <C Ai is no longer fulfilled. Then the nucleation probability per atom becomes of order 
unity, and the nucleation rate on top of the island is w « FA. Repeating the above analysis for this case, one finds 
that 



/ = 1 - exp[-6» 2 /2iV] = 1 - exp[-(L/L c 



(25) 



with L x = (2/a 2 N) 1 / 4 . The critical island size is of the order of. the square root of the island spacing, independent 
of v' and i* . This corresponds to regime I of Rottler and MaassO. 



B. Conditions for layer- by-layer growth and surfactant action 

According to classical nucleation thcoryQ the first layer island density depends on flux and in-layer mobility 
according to the scaling law 

N~ {F/ V f'^ +2 \ (26) 



Provided the in-layer mobility on the island is the same as on the terrace, for i* = 1 the LAM expression (|24|) can 
therefore be written as 



^1/21^/1/7 
pi/21 

while the corresponding expression of the TDT theory! reads 



L - ^ F A/2i (LAM), (27) 
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l *~^fW* ( tdt )- ( 28 ) 

Surprisingly, a decrease of the in-layer mobility v is seen to (weakly) decrease the critical island size according to the 
LAM, but increases it according to the TDT theory. 

A decrease of the in-layer mobility is believed to play an important role in the ability j^niLcertain adsorbates (such 
as Sb on Ag(lll)) to act as surfactants, in the sense of promoting layer-by-layer growthQcj'tll. The transition from 
three-dimensional to layer-by-layer growth occuraJ when the critical island size becomes comparable to the distance 
TV -1 / 2 between first layer islands. Equating (p7|) and ( ^8|) to {v/ F) 1 ^ one finds in both cases L c ~ £es ~ v / v> : which 
shows that the condition for the onset of layer-by-layer growth does not depend on the details of the second layer 
nucleation mechanism, at least as long as ^es is the only additional length scale in the problem; in the case i* = 2 
the dimensionless number t^v introduces another scale, and the situation becomes more complicated. 



C. Adatom density at second layer nucleation 

It was mentioned in Section | that, in addition to the rate equation approagbBi, the problem of second layer 
nucleation has been treated using the concept of a critical density for nucleationEm. To see whether such a concept is 
meaningful in the present context, we can estimate the mean adatom density n c on top of a first layer island at the 
time of second layer nucleation. Using (|J) and (Q) we obtain n c = FaL c /v', and inserting the expression (|24|) yields 

n c ~ (F/*/) 6/7 A- 1/7 ~ L- 6 N-\ (29) 

This is to be [Compared to the adatom density n su b on the substrate at the time of first layer nucleation, which is of 
the order offl (Ffiy)N- 1 ~ N 2 for i* = 1. The comparison shows that n c 3> n su b provided L c is small comp ared to 



the mean island spacing, which is equivalent to the condition of strong step edge barriers (see Section [II B). Thus 
in the regime of interest, the assumption n c ss ?i su b is not satisfied, and consequently the estimates of L c obtained in 
Refs. 6 and 7 are not quantitatively accurate. 



D. Nucleation on predeposited islands: The case of Ag(lll) 

A detailed experimental study of second layer nucleation on Ag(lll) was performed by Bromann et al.i. They 
prepared arrays of approximately circular islands of uniform initial radius ro through Ostwald ripening at high 
temperatures. Given the initial coverage 9q, the island density is then N = Qo/irr 2 . Subsequently a second dose of 
coverage A9 was deposited, and the fraction / of islands with second layer nuclei was measured as a function of ro . 
Independent evidence shows that the critical nucleus size is i* n = 1 under the experimental growth conditions. Provided 
further nucleation during the second dose can be neglectedu, the island radius increases with coverage according to 
r = roy/l + 9 /do. Inserting this into (^|) and (^) and using that 9q — AO in the experiment, we find 

/ = l-cxp[-(r /f c ) 5 ] (30) 

with 

. _ / 7 v>_ 

? ' c ~ U 2 (2 7/2 -l) ' FA9 

The step edge barrier can be extracted directly by comparing the critical radii at two different temperatures T\ and 
T 2 , since according to ( fjl| ) 

f c (Ti)/f c (T 2 ) = exp^T? 1 - T^)/5k B }. (32) 

From Figure 2 of Ref. 8 we estimate that f c w 28A at 120 K and f c w 52A at 130 K. This yields E s ~ 0.42 eV, 
leading to an additional step edge barrier AE$ « 0.32 eV. The theory of TDT predicts instead that f c ~ (i/ 2 /v) 1 ^ , 
which implies 

f c {Ti)/f c {T 2 ) = exp[(2£ s - Ed)^ 1 - T^)/Ak B ] (33) 
and yields the estimate AE$ ~ 0.12 eV reported by Bromann et al.i. 
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The dramatic discrepancy between the two estimates for AEs makes it importantO to examine also the prefactor 
v'q in the expression for v' . Using (|3l| ) with Eg — 0.42 eV yields v' w 8 x 10 19 s _1 , which seems much too large to 
be physically reasonable. Smaller values of v' can be obtained if a (slight) temperature dependence of the barrier is 
permitted. Since the step edge barrier of interest here in fact constitutes an average oveiJjhe (temperature-dependent) 
step structure, such a dependence could be expected due thermal roughening of the stept3. For example, if we demand 
that the prefactor be the same as for in-layer transports, v' Q = Vq = 2 X 10 11 s _1 , then ( |3l| ) yields E$ ~ 0.21 eV 
atjT = 120 K and Es ~ 0.20 eV at T = 130 K, somewhat smaller than the estimate Es ~ 0.22 eV of Bromann et 
al.tl, while taking their value v' = 10 13 s _1 for the prefactor, yields the larger barrier Es — 0.26 eV (T = 120 K) 
and Es — 0.24 eV (T = 130 K), respectively. The effective barrier decreases with increasing temperature, which is 
consistent with the step roughening picture. Clearly a more definite statement about the values of Es and v' on the 
Ag(lll) surface requires a detailed reanalysis of the experimental data, and perhaps also additional measurements. 



E. Application to Pt(lll) islands with decorated steps 

In this section we use the LAM to analyze the dependence of v' and AEs on the CO partial pressure pco m Pt 
homoepitaxy on Pt(lll). CO exposure during growth leads to step decoration, which impairs the hopping of adatoms 
over the step edge and therefore increases AEs (for details see Ref. 13). 

Rather than following the evolution of the fraction / of islands with second layer nuclei with coverage or island size, 
we employed a method which requires only a single growth experiment within the appropriate coverage range, i.e. a 
growth experiment which creates islands partly with and partly without second layer nuclei. The practical analysis 
is performed by determination of the number of edge sites L of the largest islands without a second layer island on 
top and of the smallest islands with a second layer one on top for a statistically significant number of topographs. 
The averages of the two sets of values vary typically by less than 10 %. The mean of the two averages then gives an 
estimate of the island size L>i/ 2 defined by f(Lup) = 1/2: at this island size it is equally likely to find an island with 



or without a second layer nucleus. Using Eqs.([5 22 24|) for i* = 1, one obtains the expression 

1/1 = Tih" 4 ™^ 2 (34) 

for v' in terms of L 1 / 2 . 

In order to translate v' into the additional step edge barrier AEs, it is necessary to make assumptions on the adatom 
migration energy Ed and the attempt frequency v' for hopping over the step, edge. For the case of Pt adatoms on 
Pt(lll), Ed was found to be 0.26 eV in two independent experimental studiesEirE 2 ], and the prefactors agreed within 
a factor of two with an average of vq = 8 x 10 12 s . An estimate of u' for the clean case can be extracted from 
a Field Ion Microscopy study, in which the adatom residence time was directly measurecO. Using the relation (]?]) 
for circular islands, the prefactor obtained in Ref. 19 translates into i/J«7x 10 11 s _1 , somewhat smaller than v§. 
However, since nothing is known about the effect of CO decoration on v' , in the following we take v' to be identical 
to vq for simplicity. The growth experimental 3 ] were carried out at 400 K with a deposition rate F — 5 x 10~ 3 ML/s. 

The open circles in Fig. [l] exhibit the critical size L 1 / 2 as a function of the applied CO partial pressure during 
growth, as derived from Fig. 2 of Ref. 13. The island size for the lowest CO partial pressure (the 'clean' case) is 
only a lower bound for L>i/ 2 , since the first monolayer coalesces prior to second layer nucleation. The full symbols 
show the variation of the additional step edge barrier AEs with pco ■ The data points represented by full circles are 
determined by use of equation (|I|). Due to the applied CO partial pressure, the additional step edge barrier increases 
dramatically from 0.12 eV to 0.36 eV. For the calculation, a was adjusted between 1/18 and 1/15 to account for the 
island geometry (compare to Fig. 1 of Ref. 13). _ 

The results are compared to the step edge barrier obtained from the TDT theory!. We have adapted their expression 
for the nucleation rate on a circular island (Eq.(3) of Ref. 5) to the notation of the present paper, which yields 

,,,_( 121n2 52 ^ 1/2 y /2 12 aL 1/2 

1 \"fFNva*L 6 1/2 49 ^ J 7 v ' 1 ' 

A reasonable choiceS for the capture number in this case is 7 = 3. The resulting data is shown in Fig. [j] as full 
triangles. It is obvious that the TDT approach seriously underestimates the magnitude of the step edge barrier for 
large barriers; for example, for pco = 1-85 x 10~ 9 mbar the TDT formula yieldsL 3 ] AEs — 0.28 eV instead of 0.36 eV. 
We note that the strong barrier condition v/v' > L is valid only for the CO partial pressures in excess of 1 x 10 -10 
mbar, hence the barrier obtained using the LAM can be trusted only for these pressures. Attempts to improve the 
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estimates for small barriers using the upper bound ( jl3| ) on the nucleation rate result in negative values for v 1 ', which 
shows that the bound is too rough to be useful. 

We still need to address the consistency of our assumption that i* = 1. In fact, at 400 K a Pt-dimer may-iiot be 
consideredrto be a stable nuclcusEIJ, because Pt-dimer breaking energies on Pt(lll) obtained in experimentQO and 
calculatior£3 vary from 0.49 eV to 0.81 eV. It is likely that at 400 K a trimer is a st able nucleus. At this temperature, 
the condition r tr -C r^s -C At ■ (rt r /r) for the scaling regime III of Section |ll C| is fulfilled. For example, for the 
largest CO partial pressure and assuming a dimer breaking energy of 0.65 eV one obtains with the other values above 
Ttr = 1.4 x 10 -6 s, Tdis = 1.9 x 10~ 5 s and At ■ {rtr/r) — 2.4 x 10~ 3 s. Unfortunately the nucleation rate in regime III 
depends explicitly on the (unknown) dimer dissociation time, and therefore the expression ( |B3| ) cannot be used for a 
quantitative analysis. Since nucleation on a top terrace is more difficult for i* > 1 compared to i* — 1, and at a given 
island size a more efficient step edge barrier is necessary for nucleation to occur, the step edge barrier obtained under 
the assumption i* = 1 is only a lower bound. 

We conclude, therefore, that the step edge barriers presented in Fig. [l] for pressures pco > 10 -10 mbar are lower 
bounds on the true values. For smaller pressures this need no longer be true, because our expression (||) underestimates 
the nucleation rate. For the data point corresponding to 'clean' growth (pco = 5 x 10~ 12 mbar) the situation is further 
complicated by the fact that the measurement yields only a lower bound on Li/2- Therefore at this point a direct 
comparison with other measurements^ and calculations!!!! of the additional barrier for clean steps on Pt(lll) is not 
possible. 



IV. THE SIZE OF THE TOP TERRACE OF A WEDDING CAKE 



In this section we apply the LAM to the growth of multilayer films in the "wedding cake" regimeEEHll. Wedding 
cakes are pyramidal mounds of fixed lateral size which grow on the template of the first layer islands. Th ey constitute 
the typical growth morphology under conditions of strong step edge barriers in the sense of Section [I A . In the limit 
of infinite step edge barriers (£es — * oo) the wedding cakes have a characteristic pointed shape which can be deducedcil 
from the observationEa that the coverages of the exposed layers follow a Poisson distribution. When fgg is finite but 
large the overall shape remains unchanged, however the pointpdr-fcip is replaced by a flat top terrace, whose typical 
size L top directly reflects the strength of the step edge barrieraL§cJ. 

For an order of magnitude estimate it suffices to note that the approximate reproduction of the mound shape after 
the growth oLan additional layer requires on average one nucleation event per top terrace and deposited monolayer, 
and thereforeEj 



u(L top )/F w 1. 



(36) 



The purpose of this section is to refine this argument and to develop a simple .analytic model for the size distrib ution 
of top terraces. A detailed comparison with growth experiments on Pt(l ll)EJ wil l be c arried out in Section [VB, 
while extensions and modifications of the theory are discussed in Sections IV C and IV D 



A. Model for the dynamics of the top terrace 

Any degree of interlayer transport couples the growth of all atomic layersEL However, the observationllljll that for 
strong barriers the shape is affected only near the top of the wedding cake suggests that the main features can be 
captured in a model which decouples the growth of the top terrace from the rest of the structure. This is achieved 
most simply by assuming that the top terrace grows on a base terrace of fixed size A, which is bounded by a step with 
an infinite edge barrier. Then all atoms landing on the base terrace or on the top terrace contribute to the growth of 
the latter, and its size evolves according to 

L{t) = VYtk = V9A, (37) 

where time t and coverage 8 = Ft are counted from the event of nucleation. The interpretation of A as the size of the 
base terrace should not be taken too literally, since we will allow for L to become larger than A; rather, A sets the 
scale of the top terraces and will eventually be fixed self consistently. In principle the constraint L < A could be built 
into the growth law for L(t ), bu t we have found that this does not improve the agreement with the experimenta l data 
to be presented in Section |IVB . An example of such a modified growth law will be briefly discussed in Section IV D 



The growth law (|37|) holds until a new top terrace nucleates on the previous one, at which point L is reset to 0. 
The growth is assumed to be deterministic, which implies that fluctuations in the deposition and diffusion events on 
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the terrace are ignored. At least for sufficiently large terraces these should be negligible compared to the fluctuations 
in the times between subsequent nucleation events, which are taken into account explicitly. The probability for no 
nucleation to occur up to time t (coverage 8 — Ft) is 



P Q (6) = exp[- / dt' u(L(t'))] = exp 



_ nA fc0(fc+2)/2 



(38) 

Here the general expression (^Tj) for the nucleation rate as a function of island size and the relation fl37| ) for the time 
dependence of L have been used. The distribution P(0) of coverages between subsequent nucleation events is then 
given by P{&) = —dPg/dd. To fix the value of A we observe that, in order for the mound morphology to be reproduced 
after the growth of one layer, the mean of P(0) should equal unity. This implies the condition 

/•OO 

dJB 6P{6) = / P (6) = 1. (39) 



A= (h+lcX^tl-V* (40) 



Inserting (B8j) we find 



where 

C k = [T((k + 4)/(k + 2))}^+ 2 ^ 2 (41) 

and r denotes the T-function. For example, in the case i* — 1 we have k = 5 and Q — c?Fjv' and (^) yields 

A = 1.193 -a- 3/5 (^7F) 1/5 . (42) 

It is plausible (and can be checked by a formal argument) that the probability Q(L) of finding a top terrace of size 
L is proportional to the probability that no nucleation has occured up to coverage 8 = (L/A) 2 . Using the expression 
( |38| ) and transforming variables from 8 to L using ( j37j ) we find 

9 T 

Q(L) = (d8/dL)(8)- l P a {L 2 /K 2 ) = exp[-C* fe (L/A) fc + 2 ]. (43) 

Here the condition (8) = 1 has been enforced through the relation (fi"o|). 

The distribution increases linearly for small L and is rather abruptly cut off at L as A (see Fig. |J). Its moments 
can be expressed in terms of r-functions. In particular, the mean size of the top terrace is given by 

qwA-g r (( fc + 5 V fc + 2 )) m 

W/A "3r((fc + 4)/(fc + 2))3/2- i 44 ) 

For k — ►oo the factor on the right hand side approaches 2/3, which corresponds to a linear distribution Q(L) = 2L/A 2 
with a sharp cutoff at L = A. Already for k = 5 the ratio (L)/A w 0.692 is very close to 2/3. Another consequence of 
the skewed shape of Q(L) is that the most probable value L max is considerably larger than the mean, and is generally 
given by 

W/A = (k + 2)- 1 /( fc+2 )r((fc + 4)/(fc + 2))- 1 / 2 . (45) 
t /A w 0.80, and i max /A — > 1 for k — > oo. 



B. Multilayer growth on Pt(lll) 

The model for the dynamics of the top terrace is substantiated by its application to mound growth on Pt(lll). As 
shown in Fig. ||, after deposition of 37.1 ML Pt with a deposition rate F — 1.3 x 10~ 2 ML/s in a partial pressure 
Vco — 1-9 x 10~ 9 mbar of carbon monoxide at 440 K mounds form, which are built from terraces scattering slightly 
around the shape of an equal-sided hexagon. The differentiated representation (morphology appears illuminated from 
the left) of the scanning tunneling topographs allows to distinguish all atomic steps, in particular the steps bounding 
the top terrace and its base terrace. The applied CO partial pressure during growth results in a significant decoration 
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of step edges by CO molecules, while the equi librium coverage of the terraces is negligible at 440 K (about 7 x 10~ 3 
MLE3). As already mentioned in Section HIE , the step decoration with CO gives rise to a significant increase in the 
additional step edge barrienla. 

The distributions of the sizes L of the top terraces and A of the base terraces determined for about 150 mounds 
are presented in Figure 0. The mean base terrace size in units of edge sites is (A) = 212. The relatively narrow 
distribution of base terraces shown in Figure supports our model assumption of a constant base terrace size. The 
mean top terrace size in units of edge sites is (L) = 142. The ratio of the two averages (L) / (A) = 0.67. The most 
probable value L max is considerably larger than its mean and is found in the size interval [0.78, 0.81] (A). 

For comparison w ith th e growth model developed in the previous section a choice for the size of i* is necessary. As 
discussed in section [HE, the scaling regime III with i* — 2 most probably applies here. Consequently we perform 
the comparison of the shape of the size distribution for i* — 2, i.e. k = 6 in (J43|); note however that ( J4^) depends 
rather weakly on k, and therefore the results assuming i* = 1, k = 5 are almost identical. The experimental ratio 
(L)/(A) = 0.67 is in excellent agreement with the model prediction of (L)/A — 0.687, and the range of experimental 
values for L max agrees very well with the model prediction L max = 0.81A. The figure also shows the analytic expression 
( [43| ) for the distribution of top terrace sizes. Given the limited number of evaluated mounds, the agreement is quite 
satisfactory. 

Due to the uncertainty in the dime r bond strength, for the quantitative estimate of v' and AEs we will assume 

this will provide a lower bound on the step edge barrier. Rearrangement of 



1; as discussed in Section [HE 



formula ( |42[ ) yields 



= 0.414 • a FA . 



(46) 



Inserting the experimental deposition rate, choosing a = 1/12 for equal-sided hexagons and a triangular lattice and 
(A) = 212 as determined experime ntally , we obtain for the hopping rate over step edges v' = 1.32 x 10 6 s _1 . With 
the same assumptions as in section III E| this translates into AEs = 0.33 eV, in reasonable agreement with the step 
edge barrier obtained at 400 K and the same CO partial pressure in the analysis of second layer nucleation in Section 
HIE. The agreement indicates that in both cases the CO concentration at the step edges was close to saturation, 



and supports the consistency of our two approaches to the determination of the additional step edge barrier from 
submonolayer islands and multilayer mounds. 



C. The one-dimensional case 



The considerations of Sections [I A and IV A are easily extended to one-dimensional islands. For an island of length 
£ the mean interval time between subsequent deposition events is At = 1/F£, while the residence time is r = £/2v' . 
Therefore the nucleation rate for strong step edge barriers and i* = 1 is 



to 



(Aty 



F 2 £ 3 
2i? ' 



(47) 



The time evolution for the size of the top terrace is now £(t) = FtX = 6X, with A denoting the length of the base 
terrace, and the probability Pq{9) that no nucleation event has taken place after a coverage 9 is 



P (9)=cxp[-(FX 3 /8v')9 4 }. 
As before, A is determined by the condition (|39|). This implies 

A = (r(l/4) 4 / 3 /2 5/3 )(^/i ? ) 1/3 « 1-75 (iZ/F) 1 / 3 . 



1/3 



(48) 



(49) 



in 1 + 1 



To summarize, for i* = 1 the size of the top terrace scales with the Ehrlich-Schwoebel length as £ ES 

— 1/5 

dimensions, and as £ ES in 2 + 1 dimensions. This agrees with the results found in Refs. 15 and 25 through a less 
quantitative method. 



D. Deterministic nucleation and other modifications 



To assess the importance of the fluctuations in the nucleation times, we now consider a model where nucleation occurs 
at fixed time intervals corresponding precisely to the growth of one monolayer. The base terrace size A is determined 
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through the requirement that these nucleation times should correspond to the maximum of the distribution, i.e. rather 
than fixing the mean of P(9) we fix the position of its peak. Noting that P(9) = —dPo/d6 = F~ 1 luPq, the condition 
dPj d6 — for the maximum can be written as 

F%=«>, (50) 

where the nucleation rate (|2l| ) has been expressed in terms of the coverage using (|3~7|). Solving (|5(]) and setting the 
corresponding coverage to unity yields 

A = {k/2) 1 ' k Vl- 1 / k . (51) 

The numerical coefficients in ( pj| ) and ( |5l| ) coincide within a few percent for 5 < k < 8. Thus for most purposes 
the simpler relation ( |5l"| ) can be used. Similarly for the one-dimensional case with i* = 1 the deterministic approach 
yields 

A = 6 1/3 (i//F) 1/3 » 1.82 (j//F) 1/3 , (52) 

very close to (^9|). 

An advantage of the deterministic approach is that more complicated growth laws for the top terrace size L(t) are 
easily incorporated. Consider for example a situation where, instead of fixing the size of the base terrace, the base 
terrace is allowed to grow such that the area Aq of its uncovered part remains constant. For consistency we should 
then also allow some fraction 1 — /z of the mass deposited onto the top and the base terraces to contribute to the 
growth of the latter. This implies that the area A(9) of the top terrace evolves according to dA/d9 = fi(A + Aq), 
which yields the growth law 

A(0) = A (e^-l). (53) 

The additional parameter fj, is fixed by the requirement that nucleation should occur at 9 = 1, at which point the top 
terrace becomes the base terrace and hence A(l) = Aq or /i = In 2. Using the expression (|2l| ) for the nucleation rate 
and imposing the relation ( |50| ) at 6 = 1 then fixes Aq to be Aq — a{kfx) 2 / k Vl~ 2 / k , and the size of the base terrace at 
the time of nucleation is 

L = v / AV^=(^ln2) 1/fc r!- 1/fc , (54) 
which differs from ( |5l|) by less than 10 % for 5 < k < 8. 

V. CONCLUSIONS 

The main purpose of this paper has been to establish some quantitatively accurate relations between microscopic 
parameters and large scale morphology for homoepitaxial growth in the presence of destabilizing step edge barriers, 
and to demonstrate their usefulness for the well-studied case of Pt/Pt(lll). Perhaps surprisingly, we found that a 
quantitative analysis of nucleation is simpler with strong step edge barriers than without them: Due to the homogeneity 
of the adatom occupation probability, and the relation p nuc — P2, our expression (^) for the nucleation rate is exact in 
the strong barrier regime for any island shape. It would be most desirable to obtain exact results also for inter mediate 



and weak barriers as well as for reversible aggregation, but for the reasons described in Sections [I B and II C such an 
extension is far from straightforward. 

On a conceptual level, the reasons for the failure of the rate equation ansatz (|l]) in the present problem need to 
be better understood. We ha ve id entified one regime in-the case i* = 2 in which our analysis agrees with the rate 



equation theoryQ (see Section II C); such a regime existsO also for i* > 2, but the precise conditions for the validity 
of (Q) are not clear. In this-|Cpiatext it is worth pointing out that the standard rate equation approach also fails in 
one-dimensional nucleationEjcd. This suggests that the rate equation treatment may generally be expected to be 
problematic in situations involving low-dimensional or confined geometries. 
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APPENDIX A: CALCULATION OF P N+1 



Here we compute the probability p n +i of finding n + 1 (noninteracting) adatoms simultaneously on the island. 
Denote by fj the time interval between the arrivals of atoms i and i + 1, and by Ti the residence time of atom i, with 
z = 1, n. The first atom arrives at time t = 0, the second at time t = t\ and so on. For all n atoms to still be on 
the island when the n + l'st atom arrives, the residence times have to satisfy the conditions t\ > U + t^+i + ... + t n . 
The U and n are exponential random variables with mean At and r, respectively. Thus 

n />oo n I* oo 

Pn+1 = (r At)-" H / dti e-^/ At l[ / dn e- Ti ' T = 

i=\Jo i=1 Jt i +t t+1 +...+t ri 



(At)-" J] / dt t e- u ' M [] e -(**+*i+i+-+*»)/r = (At)-" JJ / 

I— 1 u 1 1 u 

n 

= n(l + *At/r)- 1 . 



exp[-ti(l/Af + i/r)] 



(Al) 



For r <C At this reduces to 



1 / T \™ 



(A2) 



APPENDIX B: CRITICAL ISLAND SIZE FOR J* = 2 



In the different regimes described in Section II C , the critical island size for second layer nucleation is given by the 



expressions 



-1/7 



N~ 1/J (regime I) 



(Bl) 



/V\ 1/5 

L c ~ - AT-Vio ( re gi me n) 



(B2) 



/\ 1/4 



F 



(^T dis )- 1/8 7V- 1/8 , (regime III) 



and for general i* the nucleation rate ( |l9|) for regime II yields 

'(4 + 3**) i *!\ 1/(4+3r) /^y*/(4+3D _ i 



2a 2i '+ 2 



F 



(B3) 



(B4) 
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FIG. 1. Critical island size L 1 / 2 (open circles) and the resulting additional step edge barrer AEs (full symbols) as a function 
of CO partial pressure. Full circles were calculated according to the LAM formula (M), while full triangles were obtained using 
the TDT formula (||). 



FIG. 2. Size distribution of top terraces. The full line shows a histogram obtained from the evaluation of 145 mounds, and 
the thick dotted line shows the (appropriately normalized) theoretical prediction (^). The thin dashed line is included to 
illustrate the narrowness of the distribution of base terraces. 



FIG. 3. Typical STM image of mounds appearing on the surface after the deposition of 37.1 ML Pt/Pt(lll) in a partial CO 
pressure of 1.9 x 10~ 9 mbar. Eight top terraces and the corresponding base terraces can be seen. 
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